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Abstract 

We review mechanisms of low-temperature electronic transport through a quan- 
tum dot weakly coupled to two conducting leads. Transport in this case is domi- 
nated by electron-electron interaction. At temperatures moderately lower than the 
charging energy of the dot, the linear conductance is suppressed by the Coulomb 
blockade. Upon further lowering of the temperature, however, the conductance may 
start to increase again due to the Kondo effect. We concentrate on lateral quantum 
dot systems and discuss the conductance in a broad temperature range, which in- 
cludes the Kondo regime. 
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1. Introduction 

In quantum dot devices 1 1 1 a small droplet of electron liquid, or just a few elec- 
trons are confined to a finite region of space. The dot can be attached by tun- 
neling junctions to massive electrodes to allow electronic transport across the 
system. The conductance of such a device is determined by the number of elec- 
trons on the dot N, which in turn is controlled by varying the potential on the 
gate - an auxiliary electrode capacitively coupled to the dot 1 1 1. At sufficiently 
low temperatures, N is an integer at almost any gate voltage Vg. Exceptions 
are narrow intervals of Vg in which an addition of a single electron to the dot 
does not change much the electrostatic energy of the system. Such a degeneracy 
between different charge states of the dot allows for an activationless electron 
transfer through it, whereas for all other values of Vg the activation energy for 
the conductance G across the dot is finite. The resulting oscillatory dependence 
G{Vg) is the hallmark of the Coulomb blockade phenomenon 1 1 1. The contrast 
between the low- and high-conductance regions (Coulomb blockade valleys and 
peaks, respectively) gets sharper at lower temperatures. The pattern of periodic 
oscillations in the G vs. Vg dependence is observed down to the lowest attainable 
temperatures in experiments on tunneling through small metallic islands (21 • 

Conductance through quantum dots formed in semiconductor heterostructures 
exhibits a reacher behavior 1 1 1. In larger dots (in the case of GaAs heterostruc- 
tures, such dots may contain hundreds of electrons), the fluctuations of the heights 
of the Coulomb blockade peaks become apparent already at moderately low tem- 
peratures. Characteristic for mesoscopic phenomena, the heights are sensitive to 
the shape of a dot and to the magnetic flux threading it. The separation in gate 
voltage between the Coulomb blockade peaks, and the conductance in the valleys 
also exhibit mesoscopic fluctuations. However, the pattern of sharp conductance 
peaks separating the low-conductance valleys of the G{Vg) dependence persists. 
Smaller quantum dots (containing few tens of electrons in the case of GaAs) show 
yet another feature |4|: in some Coulomb blockade valleys the dependence G(T) 
is not monotonic and has a minimum at a finite temperature. This minimum is 
similar in origin 1 5 1 to the well-known non-monotonic temperature dependence 
of the resistivity of a metal containing magnetic impurities 1 6 1 - the Kondo ef- 
fect. Typically, the valleys with anomalous temperature dependence correspond 
to an odd number of electrons in the dot. In an ideal case, the low-temperature 
conductance in such a valley is of the order of conductance at peaks surround- 
ing it. Thus, at low temperatures the two adjacent peaks merge to form a broad 
maximum. 

The number of electrons on the dot is a well-defined quantity as long the 
conductances of the junctions connecting the dot to the electrodes is small com- 
pared to the conductance quantum /h. In quantum dot devices formed in semi- 
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conductor heterostxuctures the conductances of junctions can be tuned contin- 
uously. With the increase of the conductances, the periodic pattern in G{Vg) 
dependence gradually gives way to mesoscopic conductance fluctuations. Yet, 
electron-electron interaction still affects the transport through the device. A 
strongly asymmetric quantum dot device with one junction weakly conducting, 
while another completely open, provides a good example of that L2, ^ ^ . The 
differential conductance across the device in this case exhibits zero-bias anomaly 
- suppression at low bias. Clearly, Coulomb blockade is not an isolated phe- 
nomenon, but is closely related to interaction-induced anomalies of electronic 
transport and thermodynamics in higher dimensions Jsj- 

The emphasis of these lectures is on the Kondo effect in quantum dots. We 
will concentrate on the so-called lateral quantum dot systems [1 4Ji formed by 
gate depletion of a two-dimensional electron gas at the interface between two 
semiconductors. These devices offer the highest degree of tunability, yet allow 
for relatively simple theoretical treatment. At the same time, many of the results 
presented below are directly applicable to other systems as well, including ver- 
tical quantum dots flO'll '121, Coulomb-blockaded carbon nanotubes ri2' ll3l . 
single-molecule transistors 1 14), and magnetic atoms on metallic surfaces 1 15l . 

Kondo effect emerges at relatively low temperature, and we will follow the 
evolution of the conductance upon the reduction of temperature. On the way to 
Kondo effect, we encounter also the phenomena of Coulomb blockade and of 
mesoscopic conductance fluctuations. 

2. Model of a lateral quantum dot system 

The Hamiltonian of interacting electrons confined to a quantum dot has the fol- 
lowing general form, 

i?dot = mi ^rjAsdjs + \Y1Y1 ^ijkAsdls'dks'dis- (2.1) 

s ij ss' ijkl 

Here an operator c?|^ creates an electron with spin s in the orbital state 4>i{r) (the 
wave functions are normalized according to J dr(p*{r)(pj{r) = 6ij); h^^ = h*^ 
is an Hermitian matrix describing the single-particle part of the Hamiltonian. 
The matrix elements hijki depend on the potential U{i — r') of electron-electron 
interaction, 

h^Jkl = j drdr'4>*{r)4>*{r')U{r - r')cp^{r')cpi (r). (2.2) 

The Hamiltonian (12. 1> can be simplified further provided that the quasiparticle 
spectrum is not degenerate near the Fermi level, that the Fermi-liquid theory 
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is applicable to the description of the dot, and that the dot is in the metallic 
conduction regime. The first of these conditions is satisfied if the dot has no 
spatial symmetries, which implies also that motion of quasiparticles within the 
dot is chaotic. 

The second condition is met if the electron-electron interaction within the dot 
is not too strong, i.e. the gas parameter is small, 

Ts — (kpGo)^^ ^1, ao = Kh^/e^m* (2.3) 

Here kp is the Fermi wave vector, oq is the effective Bohr radius, k is the dielec- 
tric constant of the material, and m* is the quasiparticle effective mass. 

The third condition requires the ratio of the Thouless energy Et to the mean 
single-particle level spacing SE to be large 1161 . 

g^Er/SE-^l. (2.4) 

For a ballistic two-dimensional dot of linear size L the Thouless energy Et is of 
the order of hvp/ L, whereas the level spacing can be estimated as 

SE - hvFkp/N - ff/m*L^. (2.5) 

Here vp is the Fermi velocity and TV ~ (kpL)'^ is the number of electrons in the 
dot. Therefore, 

g'^kpL'^VN, (2.6) 

so that having a large number of electrons 1 in the dot guarantees that the 

condition Eq. J2.4t is satisfied. 

Under the conditions ( 12. 3t . i2.4\ the Random Matrix Theory (for a review see, 
1171 '18"19" '201) is a good starting point for description of non-interacting 
quasiparticles within the energy strip of the width Ep about the Fermi level [I6j. 
The matrix elements hij in Eq. ( 12. H belong to a Gaussian ensemble 1191 1201 . 
Since the matrix elements do not depend on spin, each eigenvalue e„ of the matrix 
hij represents a spin-degenerate energy level. The spacings e„+i — e„ between 
consecutive levels obey the Wigner-Dyson statistics 1 19 1; the mean level spacing 
(e„+i - e„) = 5E. 

We discuss now the second term in the Hamiltonian ( 12. H . which describes 
electron-electron interaction. It turns out 1211 1221 l23l that the vast majority of 
the matrix elements hijki are small. Indeed, in the lowest order in l/g <C 1, the 
wave functions (pi (r ) are Gaussian random variables with zero mean, statistically 
independent of each other and of the corresponding energy levels 1241 : 

WMW) = ^-fFi\r-r'\), cj^,ir)cP^ir')^^-^Fi\r-r'\). (2.7) 
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Here A ^ L?' is the area of the dot, and the function F is given by 

F(r) ~ (exp(ifc • r))Fs • (2.8) 

where (. . .)fs stands for the averaging over the Fermi surface |fc| = kp. In two 
dimensions, the function F{r) decreases with r as F cx {kpr)^'^/'^ at kpr ^ 1, 
and saturates to ^ 1 at kpr ^ 1. 

The parameter (3 in Eq. ( 12. 7> distinguishes between the presence (/3 — 1) or 
absence {(3 = 2) of the time-reversal symmetry. The symmetry breaking is driven 
by the orbital effect of the magnetic field and is characterised by the parameter 

X = (*/$o)^/5, 

where $ is the magnetic flux threading the dot and $o = hc/e is the flux quan- 
tum, so that the limits x ^ 1 x ^ 1 correspond to, respectively, /? = 1 and 
(3 = 2. Note that in the case of a magnetic field H±^ applied perpendicular to the 
plane of the dot, the crossover (at x ^ 1) between the two regimes occurs at so 
weak field that the corresponding Zeeman energy B is negligible'. 

After averaging with the help of Eqs. ( I2.7t -( I2.8> . the matrix elements (12. 2> 
take the form 



h,,ki = {2Ec + Es/2)SuS,k + Es5,k5,i + A(2//3 - 



kl- 



We substitute this expression into Hamiltonian ( 12. H . and rearrange the sum over 
the spin indexes with the help of the identity 

^^^1^2^^1'^2 ^-^'l ^1 ^-^2^2 "^'^^l^l ^^2^2' (2-9) 

where cr = (a^, cr*', cr^) are the Pauli matrices. This results in a remarkably 
simple form 1221 1231 

77i„t = EcN^ - EsS^ + A{2//3 - l)f^f (2.10) 

of the interaction part of the Hamiltonian of the dot. Here 

^ = E4x., s = 5]dL^d„,, r = E<4i (2-11) 

ns nss' n 

are the operators of the total number of electrons in the dot, of the dot's spin, 
and the "pair creation" operator corresponding to the interaction in the Cooper 
channel. 



'For example, in the experiments \25i the crossover takes place at ~ 10 mT. Zeeman energy 
in such a field _B ~ 2.5 mK, which is by an order of magnitude lower than the base temperature in 
the measurements. 
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Fig. 1 . Equivalent circuit for a quantum dot connected to two leads by tunnel junctions and capaci- 
tively coupled to the gate electrode. The total capacitance of the dot C = Cl + Cu + Cg . 

The first term in Eq. ( I2.10> represents the electrostatic energy. In the conven- 
tional equivalent circuit picture, see Fig.Q] the charging energy Ec is related to 
the total capacitance C of the dot, Ec = e^/2C. For a mesoscopic {kpL ^ 1) 
conductor, the charging energy is large compared to the mean level spacing 5E. 
Indeed, using the estimates C ^ kL and Eqs. ( 12. 3> and (12. 5> . we find 

Ec/SE r-. L/ao ^ rs^N. (2.12) 

Except an exotic case of an extremely weak interaction, this ratio is large for 
N ^ 1; for the smallest quantum dots formed in GaAs heterostructures, Ec / SE ~ 
10 E|. Note that Eqs. STAi . ( IZ6l l. and ( inH imply that 

Et/Ec ~ 1/rs > 1, 

which justifies the use of RMT for the description of single-particle states with 
energies |e„| < Ec, relevant for Coulomb blockade. 

The second term in Eq. ( I2.10> describes the intra-dot exchange interaction, 
with the exchange energy Es given by 

Es = j drdr'U{r -r')F^{\r -r'\) (2.13) 

In the case of a long-range interaction the potential U here should properly ac- 
count for the screening 1231 . For rj <C 1 the exchange energy can be esti- 
mated with logarithmic accuracy by substituting U{r) — {e^ / Kr)6{aQ — r) into 
Eq. ( I2.13t (here we took into account that the screening length in two dimensions 
coincides with the Bohr radius oq), which yields 

Es ^rs\n{l/rs)5E <^5E. (2.14) 

The estimate Eq. ( I2.14> is valid only for Ts <C 1. However, the ratio Es/5E 
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remains small for experimentally relevant^ value '--^ 1 as long as the Stoner 
criterion for the absence of itinerant magnetism 1 26 1 is satisfied. This guarantees 
the absence of a macroscopic (proportional to N) magnetization of a dot in the 
ground state 1221 . 

The third term in Eq. ( I2.10> . representing interaction in the Cooper channel, 
is renormalized by higher-order corrections arising due to virtual transitions to 
states outside the energy strip of the width Et about the Fermi level. For attrac- 
tive interaction (A < 0) the renormalization enhances the interaction, eventually 
leading to the superconducting instability and formation of a gap Aa in the elec- 
tronic spectrum. Properties of very small (Aa ^ SE) superconducting grains 
are reviewed in, e.g., |27 1; for properties of larger grains (Aa ^ Ec) see |28l. 
Here we concentrate on the repulsive interaction (A > 0), in which case A is 
very small, 

6E SE 
Hep/ET) ^ InN'^ ' 

This estimate accounts for the logarithmic renormalization of A when the high- 
energy cutoff is reduced from the Fermi energy e p down to the Thouless energy 
Et |23|. In addition, if the time-reversal symmetry is lifted {f3 — 2) then the 
third term in Eq. (j^TOJ is zero to start with. Accordingly, hereinafter we neglect 
this term altogether by setting A = 0. 

Obviously, the interaction part of the Hamiltonian, Eq. (I2.10> . is invariant 
with respect to a change of the basis of single -particle states 4>i{r). Picking up 
the basis in which the first term in Eq. i2.l\ is diagonal, we arrive at the universal 
Hamiltonian II22II23I. 



Hdot=Y.^ndLdns+Ec{N-Nof -EsS\ (2.15) 

ns 

We included in Eq. ( I2.15> the effect of the capacitive coupling to the gate elec- 
trode: the dimensionless parameter A'o is proportional to the gate voltage. 

No = CgVg/e, 

where Cg is the capacitance between the dot and the gate, see Fig.^] The relative 
magnitude of fully off-diagonal interaction terms in Eq. ( 12. U (corresponding to 
hijki with all four indices different), not included in Eq. J2.15I I. is of the order of 
1/(7 N~^^^ ^ 1. Partially diagonal terms (two out of four indices coincide) 
are larger, of the order of \f\Jg ^ A^^^/"*, but still are assumed to be negligible 
as Af > 1. 



For GaAs (m* 0.07me, k 13) the effective Bohr radius ao ~ 10 nm, whereas a typical 
density of the two-dimensional electron gas, n ~ 10^^ cm~^ 0, coiTesponds to kp = sJI-kti ~ 
10^ cmT^ . This gives kpao ~ 1. 
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As discussed above, in this limit the energy scales involved in Eq. (I2.15> form 

a well-defined hierarchy 



If all the single-particle energy levels e„ were equidistant, then the spin S of 
an even-iV state would be zero, while an odd-iV state would have S = 1/2. 
However, the level spacings are random. If the spacing between the highest oc- 
cupied level and the lowest unoccupied one is accidentally small, than the gain 
in the exchange energy, associated with the formation of a higher-spin state, may 
be sufficient to overcome the loss of the kinetic energy (cf. the Hund's rule in 
quantum mechanics). For Es ^ 5E such deviations from the simple even-odd 
periodicity are rare |22 29. JOJ. This is why the last term in Eq. ( I2.15> is often 
neglected. Eq. ( I2.15> then reduces to the Hamiltonian of the Constant Interaction 
Model, widely used in the analysis of experimental data |T1. Finally, it should be 
emphasized that SU(2)-invariant Hamiltonian ( I2.15> describes a dot in the ab- 
sence of the spin-orbit interaction, which would destroy this symmetry. 

Electron transport through the dot occurs via two dot-lead junctions. In a 
typical geometry, the confining potential forming a lateral quantum dot varies 
smoothly on the scale of the Fermi wavelength. Hence, the point contacts con- 
necting the quantum dot to the leads act essentially as electronic waveguides. 
Potentials on the gates control the waveguide width, and, therefore, the number 
of electronic modes the waveguide support; by making the waveguide narrower 
one pinches the propagating modes off one-by-one. Each such mode contributes 
2e^/ft. to the conductance of a contact. The Coulomb blockade develops when 
the conductances of the contacts are small compared to 2e^//i, i.e. when the 
very last propagating mode approaches its pinch-off |3r'"321. Accordingly, in the 
Coulomb blockade regime each dot-lead junction in a lateral quantum dot system 
supports only a single electronic mode |33 1. 

As discussed below, for Ec ^ SE the characteristic energy scale relevant 
to the Kondo effect, the Kondo temperature Tk, is small compared to the mean 
level spacing: Tk ^ 5E. This separation of the energy scales allows us to 
simplify the problem even further by assuming that the conductances of the dot- 
lead junctions are small. This assumption will not affect the properties of the 
system in the Kondo regime. At the same time, it justifies the use of the tunneling 
Hamiltonian for description of the coupling between the dot and the leads. The 
microscopic Hamiltonian of the system can then be written as a sum of three 
distinct terms. 



Es <^SE <^ Ec- 



(2.16) 




(2.17) 



which describe free electrons in the leads, isolated quantum dot, and tunneling 
between the dot and the leads, respectively. The second term in Eq. (I2.17> . the 
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Hamiltonian of the dot i/dot, is given by Eq. ( I2.15> . We treat the leads as reser- 
voirs of free electrons with continuous spectra ^fe, characterized by constant den- 
sity of states V, same for both leads. Moreover, since the typical energies lo < Eq 
of electrons participating in transport through a quantum dot in the Coulomb 
blockade regime are small compared to the Fermi energy of the electron gas in 
the leads, the spectra can be linearized near the Fermi level, = vpk; here 
k is measured from kp- With only one electronic mode per junction taken into 
account, the first and the third terms in Eq. ( I2.17> have the form 

iJloads Xl^fc^^fe^'^afes' 6 = (2.18) 

Oiks 

i?tun„cling = ^anClks^ns + H-C. (2.19) 



cxkns 



Here tan are tunneUng matrix elements (tunneling amplitudes) "connecting" the 
state n in the dot with the state k in the lead a (a = R, L for the right/left lead). 

Tunneling leads to a broadening of discrete levels in the dot. The width Tan 
that level n acquires due to escape of an electron to lead a is given by 

ran = TTiy\tl^\ (2.20) 

Randomness of single-particle states in the dot translates into the randomness of 
the tunneling amplitudes. Indeed, the amplitudes depend on the values of the 
electron wave functions at the points of the contacts, tan (t>n{fa)- For 
It'l — ^_r| ~ kpL ^ 1 the tunneling amplitudes [and, therefore, the widths 
( I2.20H in the left and right junctions are statistically independent of each other. 
Moreover, the amplitudes to different energy levels are also uncorrected, see 
Eq. dTTt : 

Pq Fa 

^an^a'n' ~ ^aa'^mi'i ^an^a'n' ~ ~~ ^f3,l^aa'^nn' i (2-21) 



The average value = Tan of the width is related to the conductance of the 
corresponding dot-lead junction 

4e2 F 

In the regime of strong Coulomb blockade {Ga ^ /h), the widths are small 
compared to the level spacing, F^ ^ 5E, so that discrete levels in the dot are 
well defined. Note that statistical fluctuations of the widths Tan are large, and 
the corresponding distribution function is not Gaussian. Indeed, using Eqs. ( I2.20> 
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and (I2.21> it is straightforward II19II20I to show that 

e-7/2 



P(^) = 5 (7 - r<,„/r„) = { ' (2.23) 

This expression is known as Porter- Thomas distribution II34I . 



3. Thermally-activated conduction 

At high temperatures, T ^ Eq, charging energy is negHgible compared to the 
thermal energy of electrons. Therefore the conductance of the device in this 
regime Goo is not affected by charging and, independently of the gate voltage, is 
given by 

111 
7^ = 7^ + 7^- (3.1) 

Lfoo '^L Lr/f 

Dependence on A^o develops at lower temperatures, T < Eq- It turns out that the 
conductance is suppressed for all gate voltages except narrow regions (Coulomb 
blockade peaks) around half-integer values of iVo. We will demonstrate this now 
using the method of rate equations L35..36J . 



3.1. Onset of the Coulomb blockade oscillations 

We start with the regime of relatively high temperatures, 

SE-^iT-^Ec, (3.2) 

and assume that the gate voltage is tuned sufficiently close to one of the points of 
charge degeneracy, 

\No~NS\<T/Ec (3.3) 

(here Nq is a half-integer number). 

Condition ( I3.2t enables us to treat the discrete single-particle levels within the 
dot as a continuum with the density of states 1 /SE. Condition Eq. ( 13. 3t . on the 
other hand, allows us to take into account only two charge states of the dot which 
are almost degenerate in the vicinity of the Coulomb blockade peak. For iVg 
close to Nq these are the state |0) with N = A^q — 1/2 electrons on the dot, and 
the state |1) with A^ = A^q* + 1/2 electrons. According to Eqs. Ost and ( l3Jt . 
the difference of electrostatic energies of these states (the energy cost to add an 
electron to the dot) is 

E+iNo) - £^11) - S|o) = 2EciN* - A^o) < T. (3.4) 
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In addition to the constraints ( 13. 2> and (13. 3> . we assume here that the inelastic 
electron relaxation rate within the dot I/te is large compared to the escape rates 
Fq,/ h. In other words, transitions between discrete levels in the dot occur before 
the electron escapes to the leads^. Under this assumption the tunnelings across 
the two junctions can be treated independently of each other (this is known as 
sequential tunneling approximation). 

With the help of the Fermi golden rule the current from the lead a into the 
dot can be written as 

= ^ IT 5^ kL| '5(6 + eV^ - 6„ - E+) (3.5) 



h 

kns 



:{^o/(a)[l - /(en)] - ^i/(e„)[l - /(Cfc)]} 



Here Vi is the probability to find the dot in the charge states \i) {i = 0, 1), 
f{oj) — [exp(ci;/T) + 1]^^ is the Fermi function, and Va is the electric potential 
on the lead a, see Fig. In writing Eq. i3.5\ we assumed that the distribution 
functions f{^k) and /(e„) are not perturbed. This is well justified provided that 
the relaxation rate exceeds the rate ~ Goc\Vl — V/?|/e at which electrons 
pass through the dot. Replacing the summations over n and k in Eq. (13. 5> by 
integrations over the corresponding continua, and making use of Eqs. ( I2.20> and 
(I2.22> . we find 

/a = — [VoF{E+ - eV^) - ViF{eV^ - £+)] , F{u) = . (3.6) 

In the steady state, the currents across the two junctions satisfy 

I = Il= -Ir. (3.7) 



Equations (13. 6> and (13. 7> . supplemented by the obvious normalization condition 
^0+7^1 = 1. allow one to find the probabilities Vi and the current across the 
dot / in response to the applied bias V = — Vr. This yields for the linear 
conductance across the dot l35l 



G = hm dl/dV = Goo 3i,b[2g^(^^_^.)/r] ■ (3-8) 

Here iVp = corresponds to the Coulomb blockade peak. At each peak, the 
conductance equals half of its high-temperature value Goo, see Eq. ( 13. U . On the 
contrary, in the Coulomb blockade valleys (iVg ^ N^), the conductance falls off 

^Note that a finite inelastic relaxation rate requires inclusion of mechanisms beyond the 
model Eq. 12.151 . e.g., electron-phonon collisions. 
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exponentially with the decrease of temperature, and all the valleys behave exactly 
the same way. Note that the sequential tunneling approximation disregards any 
interference phenomena for electrons passing the dot. Accordingly, the result 
Eq. ( 13. 8> is insensitive to a weak magnetic field. 

3.2. Coulomb blockade peaks at low temperature 

At temperatures below the single-particle level spacing in the dot 5E, the ac- 
tivation energy for electron transport equals the difference between the ground 
state energies of the Hamiltonian (I2.15> corresponding to two subsequent (in- 
teger) eigenvalues of N . Obviously, this difference includes, in addition to the 
electrostatic contribution E+{No), see Eq. i3.4\ . also a finite (and random) level 
spacing. As a result, the distance in A'o between adjacent Coulomb blockade 
peaks is no longer 1, but contains a small fluctuating contribution of the order of 
SE/Ec- Mesoscopic fluctuations of spacings between the peaks are still subject 
of a significant disagreement between theory and experiments. We will not con- 
sider these fluctuations here (see 1 18 1 for a recent review), and discuss only the 
heights of the peaks. 

We concentrate on the temperature interval 



which extends to lower temperatures the regime considered in the previous sec- 
tion, see Eq. i3.2\ . and on the gate voltages tuned to the vicinity of the Coulomb 
blockade peak, see Eq. ( 13. 3> . Just as above, the latter condition allows us to ne- 
glect all charge states except the two with the lowest energy, |0) and Due 
to the second inequality in Eq. (13. 9> . the thermal broadening of single-particle 
energy levels in the dot can be neglected, and the states |0) and |1) coincide with 
the ground states of the Hamiltonian (I2.15> with, respectively, N = A^q — 1/2 
and N = Nq +1/2 electrons in the dot. To be definite, consider the case when 



with N being an even integer; for simplicity, we also neglect the exchange term 
in Eq. ( I2.15> . Then |0) (with even number of electrons N) is the state in which 
all single particle levels below the Fermi level {n < 0) are doubly occupied. 
This state is, obviously, non-degenerate. The state |1) differs from |0) by an 
addition of a single electron on the Fermi level n — 0. The extra electron may 
be in two possible spin states, hence |1) is doubly degenerate; we denote the 
two components of |1) by |s) with s =1, |. As discussed below, the degeneracy 
eventually gives rise to the Kondo effect. However, at T 3> [see Eq. ( 13. 9H the 
quantum coherence associated with the onset of the Kondo effect is not important. 



r„ < T < 6E, 



(3.9) 



=N + l/2 



(3.10) 
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and the rate equations approach can still be used to study the transport across the 
dot J36l. 

Applying the Fermi golden rule, we write the contribution of electrons with 
spin s to the electric current I^s from lead a to the dot as 

/aa = e ^ ^ + - ^" - ^+) {^0/(a) -Vs[l- fiik)]] 

k 

We now neglect ep as it is small compared to _E+ (thereby neglecting the meso- 
scopic fluctuations of the position of the Coulomb blockade peak) and replace 
the summation over k by an integration. This yields 

/as = ^ r„o{^o/ {E+ - eV^) - Vsf (el4 - ^+)}. (3.11) 

In the steady state the currents I^s satisfy 

Ils = -Irs=I/'2 (3.12) 

(here we took into account that both projections of spin contribute equally to the 
total electric current across the dot /). Solution of Eqs. ( 13. lU and ( I3.12> subject 
to the normalization condition Vo + + Vi = I results in |23,| 

The case of odd in Eq. ( 13.101 1 is also described by Eq. ( I3.13> after replacement 
E+{No) ^ -E+iNo). 

There are several differences between Eq. ( I3.13> and the corresponding ex- 
pression Eq. ( 13. 8> valid in the temperature range ( 13. 2> . First of all, the maximum 
of the conductance Eq. ( I3.13t occurs at the gate voltage slightly (by an amount of 
the order of T/Eq) off the degeneracy point Nq — Nq, and, more importantly, 
the shape of the peak is not symmetric about the maximum. This asymmetry is 
due to correlations in transport of electrons with opposite spins through a single 
discrete level in the dot. In the maximum, the function ( I3.13> takes value 

Note that Eqs. ( I3.13> and ( I3.14t depend on the widths Fqo of the energy level 
71 = rather then on the averages over many levels in the dot, as in Eq. ( 13. 8> . 
As already discussed in Sec. |2j the widths Fqo are related to the values of the 
electron wave functions at the position of the dot-lead contacts, and, therefore, are 



-df /duj 



1 + 



(3.13) 
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random. Accordingly, the heights Gpcak of the Coulomb blockade peaks exhibit 
strong mesoscopic fluctuations. In view of Eq. M.Til . the distribution function 
of Gpcak, see Eq. ( I3.14t . is expected to be broad and strongly non-Gaussian, 
as well as very sensitive to the magnetic flux threading the dot. This is indeed 
confirmed by calculations |37, 38 1 and agrees with experimental data 1 25 , 1391 . 
The expression for the distribution function is rather cumbersome and we will 
not reproduce it here, referring the reader to the original papers 1371 1381 and 
reviews 1 18 W/^Sl) instead. 

An order-of-magnitude estimate of the average height of the peak can be ob- 
tained by replacing Tao in Eq. ( I3.14t by Fq, see Eq. M.22\ . which yields 



Gpcak ^ Goo (3.15) 

This is by a factor ^E'/T larger than the corresponding figure Gpcak = Goo/2 for 
the temperature range ( 13. 2> . and may even approach the unitary limit e'^//i) 
at the lower end of the temperature interval j3.9> . Interestingly, breaking of time- 
reversal symmetry results in an increase of the average conductance |23|. This 
increase is analogous to negative magnetoresistance due to weak localization in 
bulk systems 140] , with the same physics involved. 



4. Activationless transport through a blockaded quantum dot 

According to the rate equations theory 1351 . at low temperatures, T ^ Ec, 
conduction through the dot is exponentially suppressed in the Coulomb block- 
ade valleys. This suppression occurs because the process of electron transport 
through the dot involves a real transition to the state in which the charge of the 
dot differs by e from the thermodynamically most probable value. The proba- 
bility of such fluctuation is proportional to cxp (— E'cjiVg — Nq\/T), which ex- 
plains the conductance suppression, see Eq. \2>.1\ . Going beyond the lowest-order 
perturbation theory in conductances of the dot-leads junctions Gq allows one to 
consider processes in which states of the dot with a "wrong" charge participate 
in the tunneling process as virtual states. The existence of these higher-order 
contributions to the tunneling conductance was envisioned already in 1968 by 
Giaever and Zeller |41 1. The first quantitative theory of this effect, however, was 
developed much later |42 1. 

The leading contributions to the activationless transport, according to 1421 . 
are provided by the processes of inelastic and elastic co-tunneling. Unlike the 
sequential tunneling, in the co-tunneling mechanism, the events of electron tun- 
neling from one of the leads into the dot, and tunneling from the dot to the other 
lead occur as a single quantum process. 
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4.1. Inelastic co-tunneling 

In the inelastic co-tunneling mechanism, an electron tunnels from a lead into one 
of the vacant single-particle levels in the dot, while it is an electron occupying 
some other level that tunnels out of the dot, see Fig. |2ja). As a result, transfer 
of charge e between the leads is accompanied by a simultaneous creation of an 
electron-hole pair in the dot. 




(b) 




Fig. 2. Examples of the co-tunneling processes. 

(a) inelastic co-tunneling: transferring of an electron between the leads leaves behind an electron-hole 
pair in the dot; (b) elastic co-tunneling; (c) elastic co-tunneling with a flip of spin. 



Here we will estimate the contribution of the inelastic co-tunneling to the con- 
ductance deep in the Coulomb blockade valley, i.e. at almost integer A^q- Con- 
sider an electron that tunnels into the dot from the lead L. If energy lo of the 
electron relative to the Fermi level is small compared to the charging energy, 
uj ^ Ec, then the energy of the virtual state involved in the co-tunneling process 
is close to Ec- The amplitude of the inelastic transition via this virtual state to 
the lead R is then given by 

An,m = (4.1) 

Ec 

The initial state of this transition has an extra electron in the single-particle state 
k in the lead L, while the final state has an extra electron in the state k' in the lead 
R and an electron-hole pair in the dot (state n is occupied, state m is empty). The 
conductance is proportional to the sum of probabilities of all such processes, 

g2 g2 2 

n,m ^ n,m 

where we made use of Eq. (I2.20> . Now we estimate how many terms contribute 
to the sum in ( 14. 2> . Given the energy of the initial state uj, the number of available 
final states can be found from the phase space argument, familiar from the calcu- 
lation of the quasiparticle lifetime in the Fermi liquid theory |43 1. For u SE 



18 



L.I. Glazman and M. Pustilnik 



the energies of the states n and m lie within a strip of the width ^ lo/SE about 
the Fermi level. The total number of the final states contributing to the sum in 
Eq. ( 14. 2> is then of the order of (iv/SE)'^. Since the typical value of u is T, the 
average value of the inelastic co-tunneling contribution to the conductance can 
be estimated as 

2 ^2 



Gin ^ I ^r-^ 1 

Using now Eq. j2.22L we find 1421 



SE J h Ef 



c 



The terms entering the sum in Eq. j4.2t are random and uncorrected. As it 
follows from Eq. (I2.23> . fluctuation of each term in the sum is of the order of its 
average value. This leads to the estimate for the fluctuation (5Gi„ = Gin — Gin 
of the inelastic co-tunneling contribution to the conductance 

T \^ /e2 ViVn\^ f 5E 



Accordingly, fluctuation of Gin is small compared to its average. 

A comparison of Eq. (14. 3t with the result of the rate equations theory (13. 8> 
shows that the inelastic co-tunneling takes over the thermally-activated hopping 
at moderately low temperatures 



T < Tin = Ec 



In 



Gl + Gr 



(4.5) 



The smallest energy of an electron-hole pair is of the order of 5E. At tempera- 
tures below this threshold the inelastic co-tunneling contribution is exponentially 
suppressed. It turns out, however, that even at much higher temperatures this 
mechanism becomes less effective than the elastic co-tunneling. 



4.2. Elastic co-tunneling 

In the process of elastic co-tunneling, transfer of charge between the leads is not 
accompanied by the creation of an electron-hole pair in the dot. In other words, 
occupation numbers of single-particle energy levels in the dot in the initial and 
final states of the co-tunneling process are exactly the same, see Fig.|2jb). Close 
to the middle of the Coulomb blockade valley (at almost integer A^'o) the average 
number of electrons on the dot, N w iVo, is also an integer. Both an addition and 
a removal of a single electron cost Ec in electrostatic energy, see Eq. ( I2.15t . The 
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amplitude of the elastic co-tunneling process in which an electron is transferee! 
from lead L to lead R can then be written as 



Aei ^ yt^^tji,^— — — — - (4.6) 



The two types of contributions to the amplitude A^i are associated with virtual 
creation of either an electron if the level n is empty (e„ > 0), or of a hole if the 
level is occupied (e„ < 0); the relative sign difference between the two types of 
contributions originates in the fermionic commutation relations. 

As discussed in Sec. |2] the tunneling amplitudes tan entering Eq. (14. 6> are 
Gaussian random variables with zero mean and variances given by Eq. ( I2.21> . It 
is then easy to see that the second moment of the amplitude Eq. (14. 6> is given by 

^ ' n 

Since for Eq ^ 5E the number of terms making significant contribution to the 
sum over n here is large, and since the sum is converging, one can replace the 
summation by an integral which yields 

Substitution of this expression into 

Gei = ^—\Ali\ (4.8) 
and making use of Eq. \2.22l gives 1421 

-r- GlGr SE 

for the average value of the elastic co-tunneling contribution to the conductance. 

This result is easily generalized to gate voltages tuned away from the middle 
of the Coulomb blockade valley. The corresponding expression reads 

GlGr SE f 1 1 \ 

^ ^JhE^ l^^W^ + A^o*-^o + lj • ^^-^^^ 

and is valid when Nq is not too close to the degeneracy points iVo = Nq and 
TVo = iVg + 1 {N^ is a half-integer number): 



i{|iVo - n;\ , \No - N* - 11} » SE/Ec 
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Comparison of Eq. ( 14. 9> with Eq. ( 14. 3t shows that the elastic co-tunneUng 
mechanism dominates the electron transport akeady at temperatures 



T<T,i^ ^EcSE , (4.11) 

which may exceed significantly the level spacing. However, as we will see shortly 
below, mesoscopic fluctuations of Gei are strong |44 |, of the order of its average 
value. Thus, although Gei is always positive, see Eq. ( I4.10> . the sample-specific 
value of Gel for a given gate voltage may vanish. 

The key to understanding the statistical properties of the elastic co-tunneUng 
contribution to the conductance is provided by the observation that there are 
many (~ Ec/&E ^ 1) terms making significant contribution to the ampli- 
tude Eq. ( 14. 6> . All these terms are random and statistically independent of each 
other. The central limit theorem then suggests that the distribution of A^i is Gaus- 
sian |23|, and, therefore, is completely characterised by the first two statistical 
moments, 

^ = ^, A ,A , = A*, A*, =5f,.A ,A* (4.12) 

eL eP e/ e£ e£ el P,i el el ' 



with A^iA^i given by Eq. ( 14. 7t . This can be proven by explicit consideration of 
higher moments. For example. 



\A%\ = 2[a^,AI,) +S\A%\. (4.13) 



The non-Gaussian correction here, S jAgJ ^ l^lil (SE/Ec), is by a factor of 
6E/Ec ^ 1 smaller than the main (Gaussian) contribution. 

It follows from Eqs. ( 14. 8> . (I4.12> . and ( I4.13> that the fluctuation of the con- 
ductance SGei = Gel — Gel Satisfies 

2 



''Gd = i;(G.,) . (4.14) 



Note that breaking of time reversal symmetry reduces the fluctuations by a factor 
of 2, similar to conductance fluctuations in bulk systems, whereas the average 
conductance (|^TOJ is not affected by the magnetic field. 

It is clear from Eq. ( I4.14> that the fluctuations of the conductance are of the 
order of the conductance itself, despite naive expectations that the large num- 
ber of the contributing states results in self-averaging. The reason is that one 
has to add amplitudes, rather than probabilities, in order to compute the con- 
ductance. Because the fluctuations of the conductance are large, its distribution 
function is not Gaussian. Given the statistics (I4.12> of the amplitudes, it is not 
quite surprising that the distribution of Gei normalized to its average coincides 
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with Porter-Thomas distribution (I2.23> . This result was obtained first in 1441 by 
a different (and more general) method. 

Finally, it is interesting to compare the elastic co-tunneling contribution to 
the conductance fluctuations Eq. ( I4.14> with that of inelastic co-tunneling, see 
Eq. ( 14 ■4> . Even though the inelastic co-tunneling is the main conduction mecha- 
nism at T > Te;, see Eq. ( 14.111 1. elastic co-tunneling dominates the fluctuations 
of the conductance throughout the Coulomb blockade regime T < Eq- 



5. Kondo regime in transport through a quantum dot 

In the above discussion of the elastic co-tunneling we made a tacit assumption 
that all single-particle levels in the dot are either empty or doubly occupied. This, 
however, is not the case when the dot has a non-zero spin in the ground state. A 
dot with odd number of electrons, for example, would necessarily have a half- 
integer spin S. In the most important case of — 1/2 the top-most occupied 
single-particle level is filled by a single electron and is spin-degenerate. This 
opens a possibility of a co-tunneling process in which a transfer of an electron 
between the leads is accompanied by a flip of electron's spin with simultaneous 
flip of the spin on the dot, see Fig.|2jc). 

The amplitude of such a process, calculated in the fourth order in tunneling 
matrix elements, diverges logarithmically when the energy lo of an incoming 
electron approaches 0. Since uj ^ T, the logarithmic singularity in the trans- 
mission amplitude translates into a dramatic enhancement of the conductance G 
across the dot at low temperatures: G may reach values as high as the quantum 
limit 2e^/h 1451 1461 . This conductance enhancement is not really a surprise. In- 
deed, in the spin-flip co-tunneling process a quantum dot with odd N behaves 
as 5 = 1/2 magnetic impurity embedded into a tunneling barrier separating 
two massive conductors 1471 . It is known 1481 since mid-60's that the presence 
of such impurities leads to zero-bias anomalies in tunneling conductance II49I . 
which are adequately explained 1.50. .51 J in the context of the Kondo effect ^6). 

5.7. Effective low-energy Hamiltonian 

At energies well below the threshold A ^ SE for intra-dot excitations the transi- 
tions within the (25* + l)-fold degenerate ground state manifold of a dot can be 
conveniently described by a spin operator S. The form of the effective Hamilto- 
nian describing the interaction of the dot with conduction electrons in the leads 
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is then dictated by SU(2) symmetry"*, 

Ctks act' 

with Saa' = J2kk'ss' ^Ifes (cTfis' /2) c„/fe/y. The sum over k in Eq. is re- 
stricted to l^fe I < A. The exchange amplitudes J^^, form 2x2 Hermitian matrix 
J. The matrix has two real eigenvalues, the exchange constants Ji and J2 (here- 
after we assume that Ji > J2). By an appropriate rotation in the R — L space the 
Hamiltonian i5.2i can then be brought into the form 

Heu = J2^k^lks^,ks + E ^7(^7 • S). (5.2) 

■yks 7 

Here the operators ip^ are certain linear combinations of the original operators 
l describing electrons in the leads, and 

S7 = E '^Ik.^'^'rk's' 

kk' ss' 

is local spin density of itinerant electrons in the "channel" 7 = 1,2. 

The symmetry alone is not sufficient to determine the exchange constants J^; 
their evaluation must rely upon a microscopic model. Here we briefly outline the 
derivation i33l 1521 1531 of Eq. \5.l\ for a generic model of a quantum dot system 
discussed in Section|2labove. For simplicity, we will assume that the gate voltage 
iVo is tuned to the middle of the Coulomb blockade valley. The tunneling ( I2.19> 
mixes the state with N ^ No electrons on the dot with states having ± 1 
electrons. The electrostatic energies of these states are high Ec), hence the 
transitions ^ ± 1 are virtual, and can be taken into account perturbatively 
in the second order in tunneling amplitudes 1 54 1 . 

For the Hamiltonian (I2.15> the occupations of single-particle energy levels are 
good quantum numbers. Therefore, the amplitude Jaa' can be written as a sum 
of partial amplitudes, 

Jact' — ^ ^ ^aa' • (5.3) 
n 

Each term in the sum here corresponds to a process during which an electron or 
a hole is created virtually on the level n in the dot, cf. Eq. ( 14. 6> . For Ga ^ e^/h 
and Es <C SE the main contribution to the sum in Eq. (15. 3> comes from singly- 
occupied energy levels in the dot. A dot with spin S has 25* such levels near the 

'*In writing Eq. 15. li we omitted the potential scattering terms associated witli tlie usual elastic 
co-tunneling. This approximation is well justified when the conductances of the dot-lead junctions 
are small, Got <C /h, in which case G^i is also very small, see Eq. 14.71 
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Fermi level (hereafter we assign indexes n = —S,...,n — S to these levels), 
each carrying a spin S/25, and contributing 

Jaa'^^Cjo^'n, K^2/S, \n\<S (5.4) 
to the exchange amplitude in Eq. (15. 1> . This yields 

Jaa' ~ ^ ^ ^aa' • (5.5) 
|n|<S 

It follows from Equations (15. 3> and ( 15 .4> that 

tri=^5]A„(|<i„| + |4j)- (5.6) 

By restricting the sum over n here to |n| < S*, as in Eq. i5.5\ . and taking into 
account that all A„ in Eq. (15. 4> are positive, we find Ji + J2 > 0. Similarly, from 

det J = ^ A™ A„ |I?^„ I , P,„„ = det ) (5.7) 



m.n 



and Equations (15. 4t and ( 15. 5> follows that J1J2 > for 5 > 1/2. Indeed, in this 
case the sum in Eq. (|^) contains at least one contribution with m ^ n; all such 
contributions are positive. Thus, both exchange constants J1.2 > if the dot's 
spin 5 exceeds 1/2 1.33 J. The peculiarities of the Kondo effect in quantum dots 
with large spin are discussed in Section lS^ below. 

We now turn to the most common case of S* = 1/2 on the dot j4||. The ground 
state of such dot has only one singly-occupied energy level (n = 0), so that 
det J « 0, see Eqs. i5.5\ and (15. 7> . Accordingly, one of the exchange constants 
vanishes, 

J2 « 0, (5.8) 

while the remaining one, Ji = tr J, is positive. Equation ( 15. 8> resulted, of course, 
from the approximation made in Eq. (15. 5> . For the model ( I2.15> the leading cor- 
rection to Eq. (15. 5> originates in the co-tunneling processes with an intermediate 
state containing an extra electron (or an extra hole) on one of the empty (doubly- 
occupied) levels. Such contribution arises because the spin on the level n is not 
conserved by the Hamiltonian ( I2.15> . unlike the corresponding occupation num- 
ber Straightforward calculation L52J yields the partial amplitude in the form of 
Eq. i5A\ . but with 

_ 2EcEs 
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Unless the tunneling amplitudes i^o to the only singly-occupied level in the 
dot are anomalously small, the corresponding correction 



5E ' 



to the exchange amplitude \5.5\ is small, 

Ota' 



see Eq. ( I2.16> . To obtain this estimate, we assumed that all tunneling amplitudes 
tan are of the same order of magnitude, and replaced the sum over n in Eq. ( 15 .9> 
by an integral. A similar estimate yields the leading contribution to dct J, 

or 

J2lJi^-Esl6E. (5.10) 

According to Eq. ( I5.10> . the exchange constant J2 is negative 1551 . and its 
absolute value is small compared to Ji. Hence, Eq. ( 15. 8t is indeed an excellent 
approximation for large chaotic dots with spin S = 1/2 as long as the intra-dot 
exchange interaction remains weak,^ i.e. for Es ^ 5E. Note that corrections 
to the universal Hamiltonian (I2.15> also result in finite values of both exchange 
constants, | J2I ~ JiN^^^^, and become important for small dots with N < 
10 |46 1. Although this may significantly affect the conductance across the system 
in the weak coupling regime T > Tk, it does not lead to qualitative changes 
in the results for S ~ 1/2 on the dot, as the channel with smaller exchange 
constant decouples at low energies II56I . see also Section l5^ below. With this 
caveat, we adopt the approximation ( 15. 8> in our description of the Kondo effect 
in quantum dots with spin S — 1/2. Accordingly, the effective Hamiltonian of 
the system i5.2\ assumes the "block-diagonal" form 

H,ff = Hi+H2 (5.11) 
Hi=Y.^k^lks^iks + A^i-S) (5.12) 

ks 
ks 



^Equation <5.8t holds identically for the Anderson impurity model |51] frequently employed to 
study transport through quantum dots. In that model a quantum dot is described by a single en- 
ergy level, which formally corresponds to the infinite level spacing Umit 5E ^ 00 of the Hamilto- 
nian 12.151. 
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with J = tr J > 0. 

To get an idea about the physics of the Kondo model (see [57] for recent 
reviews), let us first replace the fermion field operator Si in Eq. (I5.12> by a single- 
particle spin-1/2 operator Si. The ground state of the resulting Hamiltonian of 
two spins 

H = J(Si • S) 

is obviously a singlet. The excited state (a triplet) is separated from the ground 
state by the energy gap Ji . This separation can be interpreted as the binding en- 
ergy of the singlet. Unlike Si in this simple example, the operator Si in Eq. ( I5.12> 
is merely a spin density of the conduction electrons at the site of the "magnetic 
impurity". Because conduction electrons are freely moving in space, it is hard for 
the impurity to "capture" an electron and form a singlet. Yet, even a weak local 
exchange interaction suffices to form a singlet ground state 1 58. .59J . However, 
the characteristic energy (an analogue of the binding energy) for this singlet is 
given not by the exchange constant J, but by the so-called Kondo temperature 

Ta' - Aexp(-l/;/J). (5.14) 



Using A ^ 6E and Equations (15. 6> and (I2.22> . one obtains from Eq. ( I5.14t the 

estimate 

'SE\ 1 e^/h Ec 



Since Gq ^ /h and Ec ^ 5E, the r.h.s. of Eq. ( I5.15> is a product of two 
large parameters. Therefore, the Kondo temperature Tk is small compared to the 
mean level spacing, 

Tk < 5E. (5.16) 

It is this separation of the energy scales that justifies the use of the effective 
low-energy Hamiltonian ( 15. H . \5.2l for the description of the Kondo effect in 
a quantum dot system. The inequality ( 15.161 1 remains valid even if the conduc- 
tances of the dot-leads junctions Gq are of the order of 2e^ /h. However, in this 
case the estimate ( I5.15t is no longer applicable f60l. 

In our model, see Equations ( I5.1U -( I5!T31 . one of the channels {'ip2) of con- 
duction electrons completely decouples from the dot, while the ^/'i -particles are 
described by the standard single-channel antiferromagnetic Kondo model 161 1571 . 
Therefore, the thermodynamic properties of a quantum dot in the Kondo regime 
are identical to those of the conventional Kondo problem for a single magnetic 
impurity in a bulk metal; thermodynamics of the latter model is fully studied I 161I . 
However, all the experiments addressing the Kondo effect in quantum dots test 
their transport properties rather then thermodynamics. The electron current oper- 
ator is not diagonal in the (^i, 7/12) representation, and the contributions of these 
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two sub-systems to the conductance are not additive. Below we relate the linear 
conductance and, in some special case, the non-linear differential conductance as 
well, to the t-matrix of the conventional Kondo problem. 

5.2. Linear response 

The linear conductance can be calculated from the Kubo formula 

G=lim-^/ dte''^*([/(i),/(0)]), (5.17) 

where the current operator / is given by 

^"=li(^^"^^)' ^"^E^afesC... (5.18) 

hs 

Here is the operator of the total number of electrons in the lead a. Evaluation 
of the linear conductance proceeds similarly to the calculation of the impurity 
contribution to the resistivity of dilute magnetic alloys (see, e.g., 1621 ). In order 
to take the full advantage of the decomposition dS.lU - dsTst . we rewrite / in 
terms of the operators V'i,2- These operators are related to the original operators 
Cfl,L representing the electrons in the right and left leads via 

{ ^''^ ) = ( ]( ''''^ ). (5.19) 

V V2ks J V -Smfc'o COSfc'o / \ CLks J 

The rotation matrix here is the same one that diagonalizes matrix J of the ex- 
change amplitudes in Eq. (15. 1> : the rotation angle 6*0 satisfies the equation tan 6^ = 
ItLo/tRol With the help of Eq. j5.19> we obtain 

Nr-Nl = cos(20o)(^i - ^2) - sin(20o) ^(V^L«V^2fcs + H.c.) (5.20) 

ks 

The current operator / entering the Kubo formula ( I5.17> is to be calculated with 
the equilibrium Hamiltonian (I5.11> - (I5.13> . Since both Ni and N2 commute with 
Hcff, the first term in Eq. j5.2m makes no contribution to /. When the sec- 
ond term in Eq. ( I5.20t is substituted into Eq. ( I5.18> and then into the Kubo for- 
mula ( I5.17> . the result, after integration by parts, can be expressed via 2-particle 
correlation functions such as (■0{(t)'(/;2(t)V'2(O)'0i(O)) (see Appendix B of |63| 
for further details about this calculation). Due to the block-diagonal structure 
of Half, see Eq. ( 15.1 1> . these correlation functions factorize into products of the 
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single-particle correlation functions describing the (free) ^^2 -particles and the (in- 
teracting) i/)i-particles. The result of the evaluation of the Kubo formula can then 
be written as 

G = GoJ ^I^h'T'^ImT.H]. (5.21) 

Here 

Go = ^ sin\2e,) = ^^lo^lt^ (5.22) 

f{Lo) is the Fermi function, and (cj) is the t-matrix for the Kondo model j5.12> . 
The t-matrix is related to the exact retarded Green function of the V^i-particles in 
the conventional way, 

GksM^) = Glioj) + GI{lo)T,{uj)GI,{uj), Gl = {uj-ik + »0)-i. 

We.xeGks,k's{ijj) is the Fourier transform of Gfes,fe's(i) = ~*^(i)({^ifes W'^L^s})' 
where (...) stands for the thermodynamic averaging with the Hamiltonian \5.\2\ . 
In writing Eq. j5.2H we took into account the conservation of the total spin 
(which implies that Gks.k's' = Sss'Gks.k's, and that the interaction in Eq. ( I5.12> 
is local (which in turn means that the t-matrix is independent of k and k'). 



5.3. Weak coupling regime: Tk ^ T ^ 5E 

When the exchange term in the Hamiltonian ( I5.12> is treated perturbatively, the 
main contribution to the t-matrix comes from the transitions of the type Ii64l 

\ks,a) ^\k's\a') . (5.23) 

Here the state \ks,<j) has an extra electron with spin s in the orbital state k 
whereas the dot is in the spin state a. By SU(2) symmetry, the amplitude of 
the transition J5.23> satisfies 

A\k's\a')'^\ks,a) = (c^s's • o-ff'ff) (5.24) 

Note that the amplitude is independent of fc, k' because the interaction is local. 
However, it may depend on oj due to retardation effects. 

The transition ( I5.23> is elastic in the sense that the number of quasiparticles 
in the final state of the transition is the same as that in the initial state (in other 
words, the transition ( I5.23t is not accompanied by the production of electron- 
hole pairs). Therefore, the imaginary part of the t-matrix can be calculated with 
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the help of the optical theorem 1651 . which yields 



2 

(T s'a' 



(5.25) 



The factor 1/2 here accounts for the probability to have spin a on the dot in 
the initial state of the transition. Substitution of the tunneling amplitude in the 
form (I5.24> into Eq. (I5.25> . and summation over the spin indexes with the help 
of the identity (12. 9> result in 

~'KvlYaTs = —v'^\A^{uj)\. (5.26) 

In the leading (first) order in J one readily obtains A'^^^ = J, independently of 
oj. However, as discovered by Kondo 0, the second-order contribution A^'^^ not 
only depends on uo, but is logarithmically divergent as 0: 

A^'^\uj)^v.f ln|A/cj|. 

Here A is the high-energy cutoff in the Hamiltonian ( I5.12> . It turns out |64| that 
similar logarithmically divergent contributions appear in all orders of perturba- 
tion theory, 

z/A("'(cj) = (i.J)"[ln|A/cj|]""\ 
resulting in a geometric series 

OO OO J 

zy^(w) = = J. J^[i^J In I A/cj I]" 



l-i/Jln|A/a;| ■ 

n— 1 n— ' ' ' 

With the help of Eq. < l5.14l t this can be written as 

vA{lo) = / . (5.27) 
\n\uj/TK\ 

Substitution of Eq. (|^^} into Eq. ( 15.261 1 and then into Eq. ( I5.21l i. and evaluation 
of the integral over oj with logarithmic accuracy yield for the conductance across 
the dot ^ 

G = Go ^'^^ , T:^Tk. (5.28) 

Equation ( I5.28> is the leading term of the asymptotic expansion in powers of 
1 / ln(T/rft-), and represents the conductance in the leading logarithmic approx- 
imation. 
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Equation Eq. (I5.28> resulted from summing up the most-diverging contribu- 
tions in all orders of perturbation theory. It is instructive to re-derive it now in 
the framework of renormalization group 1 66 1 . The idea of this approach rests on 
the observation that the electronic states that give a significant contribution to ob- 
servable quantities, such as conductance, are states within an interval of energies 
of the width uj ^ T about the Fermi level, see Eq. ( I5.21> . At temperatures of the 
order of Tk, when the Kondo effect becomes important, this interval is narrow 
compared to the width of the band D = A in which the Hamiltonian J5.12t is 
defined. 

Consider a narrow strip of energies of the width 6D <C D near the edge of 
the band. Any transition ( I5.23> between a state near the Fermi level and one of 
the states in the strip is associated with high (~ D) energy deficit, and, therefore, 
can only occur virtually. Obviously, virtual transitions via each of the states in 
the strip result in the second-order correction ,P /D to the amplitude A{uo) 
of the transition between states in the vicinity of the Fermi level. Since the strip 
contains v6D electronic states, the total correction is 1 66 1 

vJ^5D/D. 

This correction can be accounted for by modifying the exchange constant in the 
effective Hamiltonian HeU which is defined for states within a narrower energy 
band of the width D ~ SD ll66l. 

^ctf = ^^V-L^fcs + Jo-snis^ ■ S), lai <D-5D, (5.29) 

ks 

, SD 

Jd-sd = Jd + ^Jd^- (5.30) 

Here Jo is the exchange constant in the original Hamiltonian. Note that the ffcff 
has the same form as Eq. (I5.12> . This is not merely a conjecture, but can be 
shown rigorously 1 59 67 1 . 

The reduction of the bandwidth can be considered to be a result of a unitary 
transformation that decouples the states near the band edges from the rest of the 
band |68|. In principle, any such transformation should also affect the opera- 
tors that describe the observable quantities. Fortunately, this is not the case for 
the problem at hand. Indeed, the angle Oq in Eq. ( I5.19> is not modified by the 
transformation. Therefore, the current operator and the expression for the con- 
ductance ( 15. 2H retain their form. 

Successive reductions of the high-energy cutoff D by small steps 6D can be 
viewed as a continuous process during which the initial Hamiltonian ( I5.12t with 
£) = A is transformed to an effective Hamiltonian of the same form that acts 
within the band of the reduced width D ^ A. It follows from Eq. ( I5.30> that the 
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dependence of the effective exchange constant on D is described by the differen- 
tial equation 1661 1671 



dC 



i^jf,, C = ln(A/i^). (5.31) 



With the help of Eq. ( I5.14> . the solution of the RG equation ( I5.31> subject to the 
initial condition Ja = J can be cast into the form 

Iii[D/Tk) 

The renormalization described by Eq. ( I5.31> can be continued until the band- 
width D becomes of the order of the typical energy |a;| ~ T for real transitions. 
After this limit has been reached, the transition amplitude A{u!) is calculated 
in lowest (first) order of perturbation theory in the effective exchange constant 
(higher order contributions are negligible for £) ~ w ), 

vA{lj) = vJd^\^\ = .'^i . 

In I w/Tif I 

Using now Eqs. (I5.26> and (I5.21> . we recover Eq. ( I5.28t . 

5.4. Strong coupling regime: T <^ Tr- 

As temperature approaches Tk, the leading logarithmic approximation result 
Eq. ( I5.28t diverges. This divergence signals the failure of the approximation. 
Indeed, we are considering a model with single-mode junctions between the dot 
and the leads. The maximal possible conductance in this case is 2e'^/h. To obtain 
a more precise bound, we discuss in this section the conductance in the strong 
coupling regime T <^ Tk- 

We start with the zero-temperature limit T = 0. As discussed above, the 
ground state of the Kondo model (I5.12> is a singlet |58 1, and, accordingly, is not 
degenerate. Therefore, the t-matrix of the conduction electrons interacting with 
the localized spin is completely characterized by the scattering phase shifts Sg 
for electrons with spin s at the Fermi level. The t-matrix is then given by the 
standard scattering theory expression |65J 

-^i/T,(0) = ^(§,-l), Ss=e^'^% (5.33) 

where §s is the scattering matrix for electrons with spin s, which for a single 
channel case reduces to its eigenvalue. Substitution of Eq. (I5.33> into Eq. (I5.21> 
yields 

G(0) ^Goij^sin^J, (5.34) 
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for the conductance, see Eq. ( 15. 2U . The phase shifts in Eqs. (I5.33> . (I5.34> are 
obviously defined only mod tt (that is, 5s and 5s + are equivalent). This ambi- 
guity can be removed if we set to zero the values of the phase shifts at J = in 
Eq. KTlS . 

In order to find the two phase shifts 5s, we need two independent relations. 
The first one follows from the invariance of the Kondo Hamiltonian ( I5.12> under 
the particle-hole transformation V^^,^ — > sV'lfc -s (here s = ±1 for spin-up/down 
electrons). The particle-hole symmetry impUes the relation for the t-matrix 

Ts{u) - -Tl^{-Lo), (5.35) 

valid at all uj and T. In view of Eq. (I5.33> . it translates into the relation for the 
phase shifts at the Fermi level (cj = 0) |69 1, 

5^+5i^ 0. (5.36) 

The second relation follows from the requirement that the ground state of the 
Hamiltonian (I5.12> is a singlet |69 1. In the absence of exchange (J = 0) and at 
T = 0, an infinitesimally weak (with Zeeman energy B — > -|-0) magnetic field 

Hb = BS^ (5.37) 

would polarize the dot's spin. Since free electron gas has zero spin in the ground 
state, the total spin in a very large but finite region of space V surrounding the 
dot coincides with the spin of the dot, {S'^)j=q — —1/2. If the exchange with 
electron gas is now turned on, J > 0, a very weak field will not prevent the 
formation of a singlet ground state. In this state, the total spin within V is zero. 
Such change of the spin is possible if the numbers of spin-up and spin-down 
electrons within V have changed to compensate for the dot's spin: SN^ — SN^ = 
1. By the Friedel sum rule, 6Ns are related to the scattering phase shifts at the 
Fermi level, 5Ns — Ss/tt, which gives 

(5| - (5| = vr. (5.38) 

Combining Eqs. (I5.36> and (I5.38> . we find 

Ss = s^- (5.39) 

Equation then yields for zero-temperature and zero-field conductance across 
the dot 1^ 

G(0) = Go. (5.40) 

Thus, the grow of the conductance with lowering the temperature is limited only 
by the value of Gq. This value, see Eq. ( I5.22> . depends only on the ratio of the 
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tunneling amplitudes |iLo/ifloh if l^iol = V-na]^ then the conductance at T = 
will reach the maximal value G — 2e^//i allowed by quantum mechanics |45|. 

As explained above, screening of the dot's spin by itinerant electrons amounts 
to ±7r/2 phase shifts for electrons at the Fermi level [lo — 0). However, the 
phase shifts at a finite u) deviate from the unitary limit value. Such deviation is to 
be expected, as the Kondo "resonance" has a finite width ~ Tk- There is also an 
inelastic component of scattering appearing at finite lo. Virtual transitions to ex- 
cited states, corresponding to a broken up Kondo singlet, induce a local repulsive 
interaction between itinerant electrons 1691 . which is the cause of inelastic scat- 
tering. This interaction enters the fixed-point Hamiltonian, applicable when the 
relevant energies are small compared to the Kondo temperature Tk. The fixed 
point Hamiltonian takes a relatively simple form 1 69 1 when written in the basis 
of electronic states that incorporate an extra ±7r/2 phase shift [see Eq. ( I5.39l l1 
compared to the original basis of Eq. \5.\2\ . In the new basis, 

i?fixed point = Y^iwlsVks " X] vLVk's + ^^2^^ ^TPi- (5-41) 

Here ps = J2kk' '■Vks'^k's ■ (^^^ colons denote the normal ordering). The form 
of the last two terms in the rh.s. of Eq. ( I5.41> is dictated by the particle-hole 
symmetry 1691 . The ratio of their coefficients is fixed by the physical requirement 
that the Kondo singularity is tied to the Fermi level |69| (i.e. that the phase shifts 
depend only on the distance u; to the Fermi level), while the overall coefficient 
1/Tk can be viewed as a precise definition of the Kondo temperature^. Note that 
by construction, the t-matrix for (^—particles T vanishes for oj = T = 0. 

The second term in the r.h.s. of Eq. ( I5.41> describes a purely elastic scattering, 
and yields a small w— dependent phase shift 1691 



6{cj) = ^ . (5.42) 

The last term in Eq. ( I5.41> gives rise to an inelastic contribution to T. Evaluation 
of this contribution in the second order of perturbation theory (see 1.62. .69J for 
details) results in 

-wr„H=*— -2 . (5.43) 

In order to obtain the t-matrix in the original basis Ts{uj), we note that the 
elastic scattering phase shifts here are obtained by simply adding ±7r/2 to 5, i.e. 

Ss{Lo)^s^+d{Lu). (5.44) 



*With this convention the RG expression Eq. 15.141 is regarded as an estimate of with loga- 
rithmic accuracy. 
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The relation between Ts{ijj) and T{uj), accounting for small inelastic term Tin, 
reads 



1 



(5.45) 



Note that the inelastic contribution enters this expression with an extra factor of 
exp(2i(5s). The appearance of this factor is a direct consequence of the unitarity 
of the scattering matrix in the presence of both elastic and inelastic scattering 
channels, see [69 1 . Taking into account Eqs. ( I5.42L ( I5.43> , and ( I5.44> . we find 



-■KvlvaTsiio) = 1 



2rp2 



Substitution of this expression into Eq. ( 15. 2H then yields 



G = Go 



1 - {ttT/TkY 



(5.46) 



(5.47) 



Accordingly, corrections to the conductance are quadratic in temperature - a typ- 
ical Fermi liquid result |69|. The weak-coupling (T ^ Tk) and the strong- 
coupling (T ^ Tk) asymptotes of the conductance have a strikingly different 
structure. Nevertheless, since the Kondo effect is a crossover phenomenon rather 
than a phase transition |57 58 59 61 1, the dependence G{T) is a smooth and 
featureless |70| function throughout the crossover region T ~ Tk- 

Finally, note that both Eqs. J5.28> and (I5.47> have been obtained here for 
the particle-hole symmetric model ( I5.12> . This approximation is equivalent to 
neglecting the elastic co-tunneling contribution to the conductance Gei- The 
asymptotes (I5.28> . ( I5.47> remain vaUd 1331 as long as Gei /Gq <C 1. The overall 
temperature dependence of the linear conductance in the middle of the Coulomb 
blockade valley is sketched in Fig.|3l 

The fixed point Hamiltonian ( I5.41> also allows one to calculate the correc- 
tions to zero-temperature conductance due to a finite magnetic field, which enters 
Eq. (I5.41> via a term Hb — J2ksi^/'^)-^ fksfks- The field polarizes the elec- 
tron gas, and the operators ps acquire non-zero ground state expectation value 
(ps) = —vBs/2. In order to calculate the phase shifts at the Fermi energy we 
replace P]Pi in the third term in the r.h.s. ( I5.41> by J2s{P-s)Ps, and £,k,£,k' in 
the second term by their values at the Fermi level ^fe = —Bs/2. As a result, 
Eq. ( I5.41> simplifies to 



ks 



kk's 



Ps 



from which one can read off the scattering phase shifts for (^—particles. 



6, = -sB/Ti 



K- 
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(5.47) 
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\ (5-28) 
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Fig. 3. Sketch of the temperature dependence of the conductance in the middle of the Coulomb 
blockade valley with S = 1/2 on the dot. The numbers in brackets refer to the corresponding 
equations in the text. 



The phase shifts in the original basis are given by |69] 

6. 

cf. Eq. J5.44> . Equation ( l5.34> then gives 



TT ~ 

s- + 5s 



B 
TV 



G — Gf) 



1 - [B/TkY 



B « Tk 



(5.48) 



Note that the effect of a finite magnetic field on zero-temperature conductance 
is very similar to that a finite temperature has on the conductance at zero field, see 
Eq. ( I5.47> . The similarity is not limited to the strong coupling regime. Indeed, 
the counterpart of Eq. ( I5.28> reads 



Lt — Lrn — o 

\n\B/TK) 



B > Tk. 



(5.49) 



5.5. Beyond linear response 



In order to study transport through a quantum dot away from equilibrium we add 
to the effective Hamiltonian dS.lU-dSJJt a term 



eV 



(5.50) 



describing a finite voltage bias V applied between the left and right electrodes. 
Here we will evaluate the current across the dot at arbitrary V but under the 
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simplifying assumption that the dot-lead junctions are strongly asymmetric: 

Gl <C Gr. 

Under this condition the angle 6*0 in Eq. ( I5.19> is small, 6*0 ~ \tLo/tm\ <C 1. 
Expanding Eq. ( I5.20> to linear order in 6*0 we obtain 

Hv{e^) = ^(^2 - iVi) + ey^o E(V'L.V'2fc. + H.c.) (5.51) 

ks 

The first term in the rh.s. here can be interpreted as the voltage bias between 
the reservoirs of 1- and 2-particles, cf. Eq. ( I5.50> . while the second term has an 
appearance of /c-conserving tunneling with very small (proportional to 6q ^ 1) 
tunneling amplitude. 

Similar to Eq. ( 15. 5H . the current operator /, see Eq. ( I5.18> . splits naturally 
into two parts, 

/ = /o + 51, 

ks 

ks 

It turns out that 61 does not contribute to the average current in the leading (sec- 
ond) order in 9q 1471 . Indeed, in this order 

{5i) = ^eyv^ f dt' {[d{t),d{f)]) a = ^vL^2fc., 

where (. . . )o denotes thermodynamic averaging with the Hamiltonian iJcff + 
Hv{0). The thermodynamic (equilibrium) averaging is well defined here, be- 
cause for 00 — the Hamiltonian iJcff + Hy conserves separately the numbers 
of 1- and 2-particles. Taking into account that {[0{t),0{t')])o depends only on 
the difference t = t — t' , we find 

(51) = ^OyVj^ £ dr ([6(r), 0(0)])^ = 0. 

The remaining contribution / ~ (^Iq'^ corresponds to tunneling current between 
two bulk reservoirs containing 1- and 2-particles. Its evaluation follows the stan- 
dard procedure ^Tl and yields i47l 
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for the differential conductance across the dot at zero temperature. Here Go 
coincides with the small (?o-liniit of Eq. (I5.22> . Using now Eqs. (I5.26> . (I5.27> . 
and (I5.46>. we obtain 



1 dl 



= < 



3 feV 
3^2/16 



I \nWeV\/TK)' 



\eV\ « Tk 
\eV\ » Tk 



(5.53) 



Thus, a large voltage bias has qualitatively the same destructive effect on the 
Kondo physics as the temperature does. The result Eq. ( I5.53> remains vaUd as 
long as r ^ |ey| ^ 5E. If temperature exceeds the bias, T ^ eV, the 
differential conductance coincides with the linear conductance, see Eqs. (I5.28> 
and (I5.47> above. 



5.6. Splitting of the Kondo peak in a magnetic field 

According to Eq. (I5.53> . the differential conductance exhibits a peak at zero bias. 
The very appearance of such peak, however, does not indicate that we are dealing 
with the Kondo effect. Indeed, even in the absence of interactions, tunneling via 
a resonant level situated at the Fermi energy would also result in a zero-bias peak 
in dl/dV. These two situations can be distinguished by considering the evolu- 
tion of the zero-bias peak with the magnetic field. In both cases, a sufficiently 
strong field (such that the corresponding Zeeman energy B appreciably exceeds 
the peaks width) splits the zero-bias peak in two smaller ones, located at zLV*. 
In the case of conventional resonant tunneling, eV* = B/2. However, if the split 
peaks are of the Kondo origin, then eV* w B. This doubling of the distance be- 
tween the split peaks is viewed by many as a hallmark of the many-body physics 
associated with the Kondo effect |4 1. 

In order to address the splitting of the zero-bias Kondo peak we add to the 
Hamiltonian (I5.12> a term 

Hb = BS,+J2 ^^Iks^iks (5.54) 

ks 

describing the Zeeman effect of the magnetic field^. We concentrate on the most 

'in general, magnetic field affects both the orbital and spin parts of the electronic wave function. 
However, as long as spin of the dot remains equal 1/2 (the opposite case is considered in the next 
Section), the orbital effect of the field is unimportant. It modifies the exchange amplitude J and the 
bandwidth A in the Hamiltonian 15.121 . which may eventually alter the values of Go and Tfc in 
Eq. 15.531 . However, the functional form of the dependence of dl /dV on 1/ at a fixed B remains 
intact. 
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interesting case of 

Tk^B-^ a. (5.55) 

The first inequality here ensures that the zero-bias peak is split, while the second 
one allows the development of the Kondo correlations to a certain extend. For 
simplicity, we also consider here a strongly asymmetric setup, as in Section lsTsl 
which allows us to use Eq. (I5.52> relating the differential conductance to the 
t-matrix. Note, however, that the results presented below remain qualitatively 
correct even when this assumption is lifted. 

As long as the bias \eV\ is small compared to B, the differential conduc- 
tance approximately coincides with the linear one, see Eq. ( I5.49> . Similarly, at 
\eV\ ^ B the field has a negligible effect, and dl/dV is given by the second 
line in Eq. (I5.53> . The presence of the field, however, significantly affects the 
dependence of dl/dV on V at \eV\ ^ B ^ T, and it is this regime that we 
address here. 

Clearly, for B ~ \eV\ ^ Tk the system is in the weak coupling Kondo 
regime, so that the perturbative RG procedure described in Section lSJl is an ade- 
quate tool to study it. Let us suppose, for simplicity, that in the absence of mag- 
netic field electrons are filhng the lower half of a band characterized by width 
2 A and constant density of states v, see Fig. |3a). In the presence of the field, 
the energy of an itinerant electron consists of the orbital part and the Zeeman 
energy Bs/2, see the second term in Eq. ( I5.54t . Condition + Bs/2 = ep 
cuts out different strips of orbital energies for spin-up and spin-down electrons, 
see Fig.|3b). To set the stage for the RG, it is convenient to cut the initial band 
asymmetrically for spin-up and spin-down electrons. The resulting band struc- 
ture, shown in Fig.|3^c), is the same for both directions of spin. Further symmetric 
reductions of the band in the course of RG leave this band structure intact. 
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Fig. 4. The energy bands for spin-up and spin-down itinerant electrons. Tlie solid lines denote the 
boundaries of the bands, the dashed lines correspond to the Fermi energy. 

(a) Band structure without Zeeman splitting; (b) same with B ^ 0; (c) same after the initial asym- 
metric cutting of the bands (Do = A — B). 

After the initial cutting of the band, the effective Hamiltonian Hdi + Hb 
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assumes the form 

H^Y. ik^l^ks + Ji^^ ■ S) + VBS^. (5.56) 

ks 

Here k is measured from the Fermi momentum of electrons with spin s, and we 
suppressed the subscript 1 in the operators tpiks etc. The Hamihonian ( I5.56> is 
defined within a symmetric band 

lai <i?o = A-S~ A, 

see Fig.|4jc). The parameter 77 in ( I5.56> is given by 

77 = 1 - i/J/2 (5.57) 

with J being the "bare" value of the exchange constant. The correction —vJ/2 
in Eq. (I5.57> is nothing else but the Knight shift 1721 . In a finite magnetic field, 
there are more spin-down electrons than spin-up ones, see Fig. EJb). This mis- 
match results in a non-zero ground state expectation value (sf) = —vB/2 of 
the operator sf in Eq. (I5.12> . Due to the exchange interaction, this induces an 
additional effective magnetic field — vJB/2 acting on the localized spin. 

The form of the effective Hamiltonian J5.56> remains invariant under the RG 
transformation. The evolution of J and 77 with the reduction of the bandwidth D 
is governed by the equations 

^-^4, (5.58) 

where ( = \n{A/D). Equation ( I5.58> was derived in Section l53l above: the 
derivation of Eq. ( 15.591 1 proceeds in a similar way. It is clear from the solution of 
Eq. ( I5.58t . see Eq. (I5.32> . that vJi) remains small throughout the weak coupling 
regime D ^ Tff . It then follows from Eqs. ( I5.57> and ( I5.59> that 77/3 is of the 
order of 1 in this range of D. The renormalization described by Eqs. ( I5.58> and 
(I5.59> must be terminated at D ~ 2B ^ Tk, so that ijeB remains smaller than 
D. At this point 

Iii(B/Tk) 

see Eq. (I5.32> . Solution of Eq. ( I5.59> in the scaling limit^, defined as 

HB/Tk) 



_HA/Tk) 



« 1, 



In this limit tlie Knight shift in Eq. 15.571 can be neglected. 
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takes the familiar form 1571 



1 



2\n{B/TK) 



(5.61) 



Now we are ready to evaluate the differential conductance. The Hamiltonian 
Eq. ( I5.56> with J and -q given by Eqs. (I5.60> and (I5.61> is defined in a sufficiently 
wide band D ~ 25 to allow for excitations with energy of the order of Zeeman 
splitting. Since the current operator remains invariant in the course of renormal- 
ization, the differential conductance can still be calculated from Eq. J5.52> . In 
the lowest order in vJb (Born approximation) we find 



1 dl 



7rVl6 



Go dV \n\BlTK) 



l + 2e{\eV\-'qBB) 



(5.62) 



The main new feature of Eq. (I5.62> compared to ( I5.53t is a threshold behavior 
at \eV\ = tibB. The origin of this behavior is apparent: for \eV\ > t/bB an 
inelastic (spin-flip) scattering channel opens up, causing a step in dl/dV. 

In the absence of Zeeman splitting, going beyond the Born approximation in 
the renormalized exchange amplitude Eq. ( I5.32> would be meaningless, as this 
would exceed the accuracy of the logarithmic RG. This largely remains true here 
as well, except when the bias is very close to the inelastic scattering threshold. 



\eV\-rjBB <B 



The existence of the threshold makes the t-matrix singular at = t/bB; such 
logarithmically divergent contribution appears in the third order in |50|. Re- 
taining of this contribution does not violate the accuracy of our approximations, 
and we find 



1 dl 
G^dV 



^2/16 



ln'{B/TK) 



1 + 20{\eV\ - r/sB) 



In- 



B 



ln{B/TK) \eV\~T]BB 



(5.63) 



We see now that the zero-bias peak in the differential conductance is indeed split 
in two by the applied magnetic field. The split peaks are located at ±.V* with 



V* = tibB 



1 



1 



2\n{B/TK) 



B. 



(5.64) 



see Eq. ( I5.61> . 
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Because of the logarithmic divergence in Eq. ( I5.63t this resuh needs a refine- 
ment. A finite temperature T ^ Tk would of course cut the divergence |50|; this 
amounts to the replacement of | \eV\ — rjBB\ under the logarithm in Eq. ( I5.63> 

by max| | \eV\ — 77s^?|, r|. This is fully similar to the cut-off of the zero-bias 
Kondo singularity in the absence of the field. 

The r = case, however, is different, and turns out to be much simpler than 
the 5 = Kondo anomaly at zero bias. The logarithmic singularity in Eq. (I5.63> 
is brought about by the second-order in Jb contribution to the scattering ampli- 
tude. This contribution involves a transition to the excited (S'- = +1/2) state of 
the localized spin. Unlike the B = case, now such a state has a finite lifetime: 
the spin of the dot may flip, exciting (a triplet) electron-hole pair of energy risB 
in the band. Such relaxation mechanism was first considered by Korringa l iTSl . 
In our case the corresponding relaxation rate is easily found with the help of the 
Fermi Golden Rule, 

2 \u\B/Tk) 

Due to the finite Korringa relaxation rate the state responsible for the logarithmic 
singularity is broadened by which cuts off the divergence 1741 in Eq. j5.63> . 
It also smears the step-like dependence on V, which we found in the Born ap- 
proximation, see Eq. ( I5.62> . As a result, the split peaks in the V^— dependence of 
dl/dV are broad and low. The width of the peaks is 

SV ^-/K o , (5.66) 

\n\B/TK)' 

and their height is given by 

G{V*)~G{0) \n[\n{B/TK)] 



G(0) HB/Tk) 



< 1. (5.67) 



Here we used the short-hand notation (iJ/dT^ = G{y). When evaluating C?(F*), 
we replaced | \eV\ ~ tjbB^ by in the argument of the logarithm in Eq. ( I5.63> . 
The right-hand side of Eq. (|^^} is small, and the higher-order in i^Jb correc- 
tions to it are negligible. 

The problem of the field-induced splitting of the Kondo peak was revisited 
many times since the works of Appelbaum |50|, see e.g. 1.75 ilT^ 1771 . 

5. 7. Kondo effect in quantum dots with large spin 

If the dot's spin exceeds 1/2, see Refs. 1781 IT9 80l, then, as discussed in Sec- 
tion lS.ll above. both exchange constants in the effective Hamiltonian j5.2l l are 
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finite and positive. This turns out to have a dramatic effect on the dependence of 
the conductance in the Kondo regime on temperature T and on Zeeman energy 
B. Unlike the case of 5 = 1/2 on the dot, see Fig.|3] now the dependences on 
T and B are non-monotonic: initial increase of G follows by a drop when the 
temperature is lowered I33II81I at S = 0; the variation of G with _B at T = is 
similarly non-monotonic. 

The origin of this peculiar behavior is easier to understand by considering 
the i?-dependence of the zero-temperature conductance |33|. We assume that 
the magnetic field H\\ is applied in the plane of the dot. Such field leads to the 
Zeeman splitting B of the spin states of the dot, see Eq. (I5.37> . but barely affects 
the orbital motion of electrons. 

At any finite B the ground state of the system is not degenerate. Therefore, 
the linear conductance at T = can be calculated from the Landauer formula 

s 

which relates G to the amplitude of scattering E>s-,rl of an electron with spin 
s from lead L to lead R. The amplitudes Ss;qq' form 2x2 scattering matrix 
Ss. In the basis of "channels", see Eq. (15. 2> . this matrix is obviously diagonal, 
and its eigenvalues exp(2z(5^s) are related to the scattering phase shifts 5^s- The 
scattering matrix in the original (i? — L) basis is obtained from 

= C/tdiag {e^'^-'^ } 

where C/ is a matrix of a rotation by an angle 9^, see Eq. ( I5.19> . The Landauer 
formula ( I5.68> then yields 

G = Goi^sin2(,5i,-52.), Go = ^sin2(20o), (5.69) 

s 

which generalizes the single-channel expression 

Equation ( I5.69> can be further simplified for a particle-hole symmetric model 
( 15. Indeed, in this case the phase shifts satisfy 5^-\ + — 0, cf. Eq. ( I5.36> . 
which suggests a representation 

5^s — sSy. (5.70) 
Substitution into Eq. (I5.69> then results in 



G = Gosm^Si-S2). 



(5.71) 
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7r/2 








B 



Fig. 5. Dependence of the scattering phase shifts at the Fermi level on the magnetic field for 5 > 1/2 
(left panel) and 5 = 1/2 (right panel). 

If spin on the dot S exceeds 1/2, then both channels of itinerant electrons par- 
ticipate in the screening of the dot's spin II56I . Accordingly, in the limit B Q 
both phase shifts 5^ approach the unitary limit value 7r/2, see Fig.|5] However, 
the increase of the phase shifts with lowering the field is characterized by two dif- 
ferent energy scales. These scales, the Kondo temperatures Ti and T2, are related 
to the corresponding exchange constants in the effective Hamiltonian j5.2> . 



so that Ti > T2 for Ji > J2. It is then obvious from Eq. j5.71l l that the con- 
ductance across the dot is small both at weak {B ^ T2) and strong {B ^ Ti) 
fields, but may become large Gq) at intermediate fields T2 B ^ Ti, see 
Fig- HI In other words, the dependence of zero-temperature conductance on the 
magnetic field is non-mono tonic. 



Fig. 6. Sketch of the magnetic field dependence of the Kondo contribution to the linear conductance 
at zero temperature. The conductance as function of temperature exhibits a similar non-monotonic 
dependence. 





This non-monotonic dependence is in sharp contrast with the monotonic in- 
crease of the conductance with lowering the field when S ~ 1/2. Indeed, in the 
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latter case it is the channel whose coupling to the dot is the strongest that screens 
the dot's spin, while the remaining channel decouples at low energies 1561 . see 
FiglH 

The dependence of the conductance on temperature G{T) is very similar^ to 
G{B). For example, for 5 = 1 one obtains 



G/Gq — 



1 
1 



1 

T'2 



n 2 



in(r/ri) HT/T2) 



r>Ti 



(5.72) 



The conductance reaches its maximun Gmax at T ~ \/TiT2. The value of G,, 



can be found analytically for Ti'3> T2. For S = 1 the result reads |33, 



Gn 



\n\Ti/T2) 



(5.73) 



Consider now a Coulomb blockade valley with N = even electrons and spin 
5 = 1 on the dot. In a typical situation, the dot's spin in two neighboring valleys 
(with iV± 1 electrons) is 1/2. Under the conditions of applicability of the approx- 
imation Eq. ( 15. 5> . there is a single non-zero exchange constant Jjv±i for each of 
these valleys. If the Kondo temperatures Tk are the same for both valleys with 
S ~ 1/2, then Jn+i = Jn-i = Jodd- Each of the two singly-occupied energy 
levels in the valley with 5 = 1 is also singly-occupied in one of the two neigh- 
boring valleys. It then follows from Eqs. J5.4t -( ls!6l that the exchange constants 
Ji_2 for 5 = 1 satisfy 



Ji 



— 2 i-^N^ 



1 + Jn-i) — Jo 



dd- 



Since both Ji and J2 are positive, this immediately implies that J1.2 < Jodd- 
Accordingly, both Kondo temperatures Ti 2 are expected to be smaller than Tk 
in the nearby valleys with S = 1/2. 

This consideration, however, is not applicable when the dot is tuned to the 
vicinity of the singlet-triplet transition in its ground state [ IT tl 79 80 1, i.e. 
when the energy gap A between the triplet ground state and the singlet excited 
state of an isolated dot is small compared to the mean level spacing 5E. In this 
case the exchange constants in the effective Hamiltonian \5.2\ should account 
for additional renormalization that the system's parameters undergo [82J when 
the high-energy cutoff (the bandwidth of the effective Hamiltonian) D is reduced 

'Note, however, t hat (f jt/i^ (t)V>J; (0)^^ (0)) ^ {i}\(t)^^ (0))(</ ),(t)V)^(0)) at finite T. 
Therefore, unlike Eq. <5.341 . there is no simple generalization of Eq. <5.2U to the two-channel case. 
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from D ^ 5E down to £> ^ A <C 5E, see also |63|. The renormalization 
enhances the exchange constants Ji 2- If the ratio A/SE is sufficiently small, 
then the Kondo temperatures Ti 2 for S = 1 may become of the same or derim 
|80|, or even significantly exceed 11 II 1121 l79l the corresponding scale Tk for 
5=1/2. 

In GaAs-based lateral quantum dot systems the value of A can be controlled 
by a magnetic field H± applied perpendicular to the plane of the dot II79I . Be- 
cause of the smallness of the effective mass m*, even a weak field H±^ has a 
strong orbital effect. At the same time, smallness of the quasiparticle g-factor in 
GaAs ensures that the corresponding Zeeman splitting remains small 1 2). The- 
ory of the Kondo effect in lateral quantum dots in the vicinity of the singlet-triplet 
transition was developed in L83J . see also I84J. 



6. Concluding remarks 

In the simplest form of the Kondo effect considered in these notes, a quantum 
dot behaves essentially as an artificial "magnetic impurity" with spin 5, coupled 
via exchange interaction to two conducting leads. The details of the temperature 
dependence G{T) of the linear conductance across the dot depend on the dot's 
spin S. In the most common case of 5 = 1/2 the conductance in the Kondo 
regime monotonically increases with the decrease of temperature, potentially up 
to the quantum limit 2e? jh. Qualitatively (although not quantitatively), this in- 
crease can be understood from the Anderson impurity model in which the dot is 
described by a single energy level. On the contrary, when spin on the dot exceeds 
1/2, the evolution of the conductance proceeds in two stages: the conductance 
first raises, and then drops again when the temperature is lowered. 

In a typical experiment jU, one measures the dependence of the differen- 
tial conductance on temperature T, Zeeman energy B, and dc voltage bias V . 
When one of these parameters is much larger than the other two, and is also large 
compared to the Kondo temperature Tk, the differential conductance exhibits a 
logarithmic dependence 



^_dI_ 



\ max{T,B,ey} 
In — 



Tk 



(6.1) 



characteristic for the weak coupling regime of the Kondo system, see Section l531 
Consider now a zero-temperature transport through a quantum dot with 5 = 1/2 
in the presence of a strong field B ^ Tk- In accordance with Eq. ( 16. H . the 
differential conductance is small compared to Gq both for eV <C B and for 
eV 3> B. However, the calculation in the third order of perturbation theory 
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in the exchange constant yields a contribution that diverges logarithmically at 
eV « ±5 150]. The divergence is reminiscent of the Kondo zero-bias anomaly. 
However, the full development of resonance is inhibited by a finite lifetime of 
the excited spin state of the dot |74|. As a result, the peak in the differential 
conductance at eV = ztpiB is broader and lower then the corresponding peak at 
zero bias in the absence of the field, see Section lsTSI 

One encounters similar effects in studies of the influence of a weak ac signal of 
frequency fl > Tk applied to the gate electrode 1 85 1 on transport across the dot. 
In close analogy with the usual photon-assisted tunneling 1. 86 J . such perturbation 
is expected to result in the formation of satellites 1871 at eV = nMl (here n 
is an integer) to the zero-bias peak in the differential conductance. Again, the 
resonances responsible for the formation of the satellite peaks are limited by the 
finite lifetime effects II88I . 

The spin degeneracy is not the only possible source of the Kondo effect in 
quantum dots. Consider, for example, a large dot connected by a single-mode 
junction to a conducting lead and tuned to the vicinity of the Coulomb blockade 
peak 1311 . If one neglects the finite level spacing in the dot, then the two almost 
degenerate charge state of the dot can be labeled by a pseudospin, while real spin 
plays the part of the channel index (^IIIH- This setup turns out to be a robust 
realization f3r'891 of the symmetric (i.e. having equal exchange constants) two- 
channel S* = 1/2 Kondo model |56 1. The model results in a peculiar temperature 
dependence of the observable quantities, which at low temperatures follow power 
laws with manifestly non-Fermi-liquid fractional values of the exponents f90T. 

It should be emphasized that in the usual geometry consisting of two leads 
attached to a Coulomb-blockaded quantum dot with S = 1/2, only the conven- 
tional Fermi-liquid behavior can be observed at low temperatures. Indeed, in this 
case the two exchange constants in the effective exchange Hamiltonian i5.2\ are 
vastly different, and their ratio is not tunable by conventional means, see the dis- 
cussion in Section lsHl above. A way around this difficulty was proposed in 1 911 . 
The key idea is to replace one of the leads in the standard configuration by a very 
large quantum dot, characterized by a level spacing 6E' and a charging energy 
Eq. At T ^ 5E', particle-hole excitations within this dot are allowed, and elec- 
trons in it participate in the screening of the smaller dot's spin. At the same time, 
as long as T ^ E'q, the number of electrons in the large dot is fixed. There- 
fore, the large dot provides for a separate screening channel which does not mix 
with that supplied by the remaining lead. In this system, the two exchange con- 
stants are controlled by the conductances of the dot-lead and dot-dot junctions. 
A strategy for tuning the device parameters to the critical point characterized by 
the two-channel Kondo physics is discussed in |92 1. 

Finally, we should mention that the description based on the universal Hamil- 
tonian ( I2.15> is not applicable to large quantum dots subjected to a quantizing 
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magnetic field H± 193II94I . Such field changes drastically the way the screening 
occurs in a confined droplet of a two-dimensional electron gas [951. The droplet 
is divided into alternating domains containing compressible and incompressible 
electron liquids. In the metal-like compressible regions, the screening is almost 
perfect. On the contrary, the incompressible regions behave very much like insu- 
lators. In the case of lateral quantum dots, a large compressible domain is formed 
near the center of the dot. The domain is surrounded by a narrow incompressible 
region separating it from another compressible ring-shaped domain formed along 
the edges of the dot |96|. This system can be viewed as two concentric capac- 
itively coupled quantum "dots" - the core dot and the edge dot 1931 1961 . When 
the leads are attached to the edge dot, the measured conductance is sensitive to 
its spin state: when the number of electrons in the edge dot is odd, the conduc- 
tance becomes large due to the Kondo effect |93|. Changing the field causes 
redistribution of electrons between the core and the edge, resulting in a striking 
checkerboard-like pattern of high- and low-conductance regions | 93 94 J . This 
behavior persists as long as the Zeeman energy remains small compared to the 
Kondo temperature. Note that compressible regions are also formed around an 
antidot - a potential hill in a two-dimensional electron gas in the quantum Hall 
regime 1971 . Both Coulomb blockade oscillations and Kondo-like behavior were 
observed in these systems 1981 . 

Kondo effect arises whenever a coupling to a Fermi gas induces transitions 
within otherwise degenerate ground state multiplet of an interacting system. Both 
coupling to a Fermi gas and interactions are naturally present in a nanoscale trans- 
port experiment. At the same time, many nanostructures can be easily tuned to 
the vicinity of a degeneracy point. This is why the Kondo effect in its various 
forms often influences the low temperature transport in meso- and nanoscale sys- 
tems. 

In these notes we reviewed the theory of the Kondo effect in transport through 
quantum dots. A Coulomb-blockaded quantum dot behaves in many aspects as 
an artificial "magnetic impurity" coupled via exchange interaction to two con- 
ducting leads. Kondo effect in transport through such "impurity" manifests itself 
in the lifting of the Coulomb blockade at low temperatures, and, therefore, can 
be unambiguously identified. Quantum dot systems not only offer a direct access 
to transport properties of an artificial impurity, but also provide one with a broad 
arsenal of tools to tweak the impurity properties, unmatched in conventional sys- 
tems. The characteristic energy scale for the intra-dot excitations is much smaller 
than the corresponding scale for natural magnetic impurities. This allows one to 
induce degeneracies in the ground state of a dot which are more exotic than just 
the spin degeneracy. This is only one out of many possible extensions of the 
simple model discussed in these notes. 
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